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ABSTRACT 

Context. The Sunyaev-Zeldovich (SZ) effect is a powerful tool for studying clusters of galaxies and cosmology. Large mm-wave 
telescopes are now routinely detecting and mapping the SZ effect in a number of clusters, measure their comptonisation parameter 
and use them as probes of the large-scale structure and evolution of the universe. 

Aims. We show that estimates of the physical parameters of clusters (optical depth, plasma temperature, peculiar velocity, non-thermal 
components etc.) obtained from ground-based multi-band SZ photometry can be significantly biased, owing to the reduced frequency 
coverage, to the degeneracy between the parameters and to the presence of a number of independent components larger than the 
number of frequencies measured. We demonstrate that low-resolution spectroscopic measurements of the SZ effect that also cover 
frequencies > 270 GHz are effective in removing the degeneracy. 

Methods. We used accurate simulations of observations with lines-of-sight through clusters of galaxies with different experimental 
configurations (4-band photometers, 6-band photometer, multi-range differential spectrometer, full coverage spectrometers) and dif- 
ferent intracluster plasma stratifications. 

Results. We find that measurements carried out with ground-based few-band photometers are biased towards high electron temper- 
atures and low optical depths, and require coverage of high frequency and/or independent complementary observations to produce 
unbiased information; a differential spectrometer that covers 4 bands with a resolution of ~ 6 GHz eliminates most if not all bias; 
full-range differential spectrometers are the ultimate resource that allows a full recovery of all parameters. 
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rays of bolometers, are now operating in excellent sites and pro- = y— — — ^-[^coth(x/2) - 4], (1) 

duce a number of detections and maps of the SZ effect, discover CMB e ' 

new clusters, and establish cluster and cosmological parameters, where I C mb is the brightness of the CMB, x = hv/kT and 

Meanwhile, the Planck space mission (Planck collaboration 

1201 lab has recently produced a shallow whole-sky survey in y - f n e cr l di 

nine cm to submm bands, from which an early catalogue of Jlos m e c 2 

massive clusters detected via the SZ effect has been extracted which is tional to the inte al of the re al 

(Planc^k collaboration |2011b|) The Early-SZ catalogue consists ±e LOS. For rich clusters the optical depth is r, = 

of 169 known clusters, plus 20 new discoveries, including ex- r ,„ , „ „, , ^ , „ „. .„ , 

, , rrn -t hi- \ ^r,^ 1 t n e o"rat ~ 0.01 and the average energy boost of CMB pho- 

ceptional members (Planck collaboration 201 lcb. JL °s „ , 

tons is AE/E ~ kT e /m e c 2 < 1%, so that a AT/T < 10~ 4 is ex- 
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CMB anisotropy, and has a characteristic spectrum, with a decre- 
ment of CMB brightness at v i 2\%GHz and an increment of 
CMB brightness at v i 2\'$>GHz. This peculiar spectrum is the 
key to distinguish the thermal SZ effect from competing contri- 
butions (see below). Relativistic corrections have also been com- 
puted for improved precision (see e.g. Itoh et al. 120001 1 and are 
used to characterize the physical parameters of the cluster atmo- 
spheres (see e.g. Colafrancesco et al. 1201 II) . 

b) A Doppler component, caused by the collective motion of 
the cluster with velocity v in the CMB restframe. The spectrum 
of this effect is the same as the spectrum of the intrinsic CMB 
anisotropy, so that the two contributions cannot be separated. 
The amplitude of the signal is 
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where Vlos is the projection of v along the LOS; so we get 
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d) the intrinsic anisotropy of the CMB (see point b) above) 
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c) A non-thermal component AI„,, caused by a non-thermal pop- 
ulation of electrons, produced by e.g. the AGNs present in the 
cluster, relativistic plasma in cluster cavities, shock accelera- 
tion. If dark matter in galaxy clusters consisted of neutralinos, 
their annihilation would produce high-energy charged particles 
as well. For a review, see e.g., Colafrancesco 120101 and refer- 
ences therein. The inverse-Compton spectrum of this component 
is very different from that of the thermal component, because the 
energy of CMB photons is boosted to frequencies much higher 
than the submm frequencies observable with CMB instruments 
(see Colafrancesco 2008 for details). This component is nor- 
mally sub-dominant with respect to the thermal component: its 
optical depth t„, is at least 50 times less than r,. The parame- 
ters characterizing the spectrum, in addition to the optical depth, 
are the spectral index of the power-law spectrum of the energy 
of the electrons a (typically around -2.7), and their minimum 
momentum p\, typically of the order of a few MeV/c. 

Additional sources of signal along the same line of sight are 
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Since components b) and d) have exactly the same spectrum, we 
will describe them in the following with the parameter AIcmb = 
A/cmbi + A/ v (or the equivalent ATcmb), which is the sum of 
the intrinsic anisotropy of the CMB and of the kinetic SZ effect 
along the line of sight. 

e) The emission of dust Aid in our Galaxy and in the galaxies 
of the cluster. This is modelled as a thermal spectrum with tem- 
perature Td ~ 2QK and a spectral index of emissivity ~ -1.5, 
or a superposition of several components with different Td- It is 
important as a contaminant at frequencies where the thermal SZ 
is positive. 

f) The free-free and synchrotron emission (A///, A/, y ) from the 
diffuse medium in our Galaxy and from the galaxies in the clus- 
ter: this component can be important at low frequencies, where 
the SZ is negative. 

According to these considerations it follows that SZ mea- 
surements promise to estimate several physical parameters of 
the cluster on the line of sight, provided there are more obser- 
vation bands than parameters to be determined, or some of the 
contributions are known to be negligible. 




400 600 
frequency (GHz) 

Fig. 1. Thermal SZ spectrum (continuous line in bottom panel), 
compared to the atmospheric transmission of a dry, cold, atmo- 
sphere (top panel, PWV=0.5mm), and to spectra of the non- 
thermal SZ effect (dotted line), of CMB anisotropy and kine- 
matic SZ (dashed line), and of dust anisotropy (dot-dashed line). 
The parameters of the different spectra are for the benchmark 
case detailed in section ^2] r, = 5 x 10~ 3 , T e = %.5keV, 
ATcmb = 22pK, t„, = 1 x 10~ 4 , a = -2.7, p x = lAMeV/c, 
A/rf(150 GHz) = 600 Jy. The frequency coverage of the dif- 
ferent experiments considered in the paper is also shown as dot- 
ted horizontal lines, labelled with the experimental configuration 
number. 



Multi-frequency measurements are therefore mandatory to 
separate the contributions of the different physical components, 
taking advantage of the characteristic spectrum of the SZ effect, 
which significantly departs from the spectra of the foreground 
and background components (see fig. [TJ. The wider and more 
detailed the frequency coverage of the observations, the more 
effective the separation of the different components. 

The recent results of Planck (Planck collaboration 1201 fbl . 
for example, have been obtained by exploiting the excellent fre- 
quency coverage of the mission and sophisticated component 
separation techniques (see e.g. Leach et al. 2008). 

In this paper we study how effective the various experimental 
configurations are in separating all the different physical compo- 
nents and in providing unbiased estimates of the cluster param- 
eters (like y, v, T e , t,, t„,, p\ ...) and of the other parameters that 
describe the signals along the same line of sight (Td, Td, I s , Iff, 

A/cMB/i ■■■)■ 

Evidently, ground-based few-band photometers cannot pro- 
vide enough information to separate all physical components. 
Observations are hampered at high frequencies ( > 200GHz) 
by atmospheric noise (see figO}: this significantly limits the 
coverage of the positive part of the thermal SZ spectrum, and 
makes the removal of parameter degeneracies much more dif- 
ficult. These instruments need external information (optical, X- 
ray, far-IR, etc.) to produce mainly measurements of r,. With the 
addition of external data, these experiments provide invaluable 
information in the current exploration phase: a large database 
of clusters is being built, and new cluster candidates have been 
discovered (see e.g. Hincks et al. 120101 Marriage et al. 2010; 
Brod win et al. EOTTl H and et al. l20TTl Sehgal etal.l20TTl Foley 
et al. I20TT1 Story et al. I20TT1 Williamson et al. l20TTb . 

In principle, future space-based spectrometers can cover the 
full range of interesting frequencies and offer much more infor- 
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mation: with these machines it should be possible to measure 
the parameters of a cluster, and use external information, when 
available, as a cross-check. Also, other important scientific tar- 
gets of these instruments are the measurement of the C + and 
CO lines in the redshift desert and beyond for a large number of 
galaxies (see e.g. de Bernardis et al. 120 101 Gong et al. 1201 11 1. 

The main purpose of this paper is to analyse quantitatively 
and compare the performance of these different configurations 
for SZ measurements. 

As samples of experimental configurations we considered a 
ground-based 4-band photometer; a 4-band photometer and a 4- 
band differential spectrometer (both suitable for balloon-borne 
observations, and similar to the configuration of the OLIMPO 
experiment - see Masi et al. 2008; Conversi et al. 1201 01 ). Planck 
HFI (Planck HFI core team 1201 ldl l. and a full-coverage differ- 
ential spectrometer suitable for a future, dedicated space mis- 
sion, like SAGACE (de Bernardis et al. 120101 ) or Millimetron 
( http ://w w w. sron.rug . nl/millimetron) . 

In section ^2] we describe and compare the considered con- 
figurations; in section |3] we describe the simulations and the 
analysis method; in section |0]we discuss the results. 



2. Experimental configurations 

All instruments considered here work at the diffraction limit, 
with AQ = A 2 , to achieve the best possible angular resolution, 
which is needed to resolve the target cluster. 

We limited our analysis to experiments using arrays of 
bolometers because of their superior mapping speed. Having 
maps of a sky patch surrounding the cluster is important to iden- 
tify foreground structures and the internal structure of the clus- 
ter, if the resolution is sufficient. This choice limits the frequency 
coverage to v i 70GHz, where bolometer arrays have a signifi- 
cant sensitivity advantage with respect to coherent detectors. 

Few-band photometers have been traditionally matched to 
the mm/submm atmospheric windows, i.e. the frequency bands 
Wl=[75-115]GHz, W2=[125-175]GHz, W3=[190-315]GHz, 
W4=[330-365]GHz, W5=[390-420]GHz. Operation in the 
higher frequency windows is significantly hampered by atmo- 
spheric noise and poor transmission, which strongly depends on 
the telescope site. In general ground-based photometers use nar- 
rower bands within the ranges above, while balloon-borne multi- 
band photometers can have wider bands and can use the higher 
frequency bands efficiently. Satellite instruments do not have 
these limitations, and the operation bands are selected based 
on the necessity to calibrate the instrument and to study the 
Galactic and extragalactic contaminating foregrounds (like the 
CO lines, and the continuum from interstellar dust, free-free and 
synchrotron emission). 

For our spectroscopy we considered a Fourier Transform 
Spectrometer (FTS) because, at variance with dispersion spec- 
trometers (see e.g. Bradford et al. 2002 ), this instrument is in- 
trinsically imaging, a crucial requirement for studying appropri- 
ately the SZ effect in cosmic structures. Moreover, at variance 
with Fabry-Perot spectrometers (see e.g. Benford et al. 120031 
) the FTS can be used in a differential configuration where a 
sky field is compared to an internal reference blackbody (as 
in the COBE-FIRAS instrument, Mather et al. 119931 or in the 
Herschel-SPIRE, Griffin et al. [20071 . and also in a differential 
configuration where two sky-fields are compared, thus rejecting 
most of the common mode signal from the instrument, the atmo- 
sphere, the foregrounds, and the CMB (see e.g. de Bernardis et 
al. 120101) . 



Experimental Configuration (ECO) is a ground-based 4- 
band photometer that measures the bands Wl, W2, and the two 
halves of W3 ([200-240]GHz and [240-3 10]GHz) in an excel- 
lent site like the Atacama desert, South Pole, or Dome-C. We 
assumed 2.7%, 2.5%, 3.5%, 5.0% average emissivity for the at- 
mosphere, respectively, in the 4 bands above (corresponding to 
a precipitable water vapor ~ 0.5 mm), and an equivalent temper- 
ature of the atmosphere of 240K. Only photon noise was con- 
sidered, assuming that the telescope is used in a differential con- 
figuration, where most of the turbulence is subtracted out as a 
common mode signal. 

Experimental Configuration 1 (EC1) is a balloon-borne 4- 
band photometer that measures bands W2, W3, W4, and W5. 
Here the background is limited by the temperature (~ 230 K) and 
emissivity of the telescope mirrors and of the cryostat window. 
Again a differential configuration was considered. 

Experimental Configuration 2 (EC2) is a balloon-borne 4- 
band differential spectrometer that measures the low-resolution 
(Av = 6 GHz) spectra in four ranges coincident with W2, W3, 
W4, and W5. Here the background is limited by the tempera- 
ture and emissivity of the telescope, of the interferometer mir- 
rors, and of the cryostat window, all at room temperature (230K). 
Note that in this (and in the following) photon-noise limited con- 
figurations, that use Fourier transform spectrometers, the error 
on each spectral bin scales as the inverse of the spectral resolu- 
tion Av (see §0. 

Experimental Configuration 3 (EC3) is an Earth-orbit (EO) 
satellite with a telescope radiatively cooled at 80K, which mea- 
sures spectra in four consecutive ranges bO = [100, 200]GHz, 
bl = [201, 350]GHz, b2 = [351, 500]GHz, b3 = [501, 600]GHz 
with 6 GHz resolution. Here the background is limited by the 
temperature and emissivity of the cold telescope and of the cryo- 
stat window (total emissivity e ~ 0.01), while the interferometer 
mirrors are kept at < 4K and do not contribute to the photon 
noise. The spectral coverage has been divided in the four ranges 
observed simultaneously by independent detector arrays in order 
to limit the radiative background on each array. 

Experimental Configuration 4 (EC4) is a 6-band photome- 
ter similar to Planck-HFI, which operates in the Lagrangian 
point L2 of the Sun-Earth system, with the telescope radia- 
tively cooled at 45K, measuring simultaneously in the bands 
bO = [83.5, 116.5]GHz, bl = [119.4, 166.6]GHz, b2 = [181.2, 
252.8]GHz, b3 = [303.6, 402.4]GHz, b4 = [460.5, 629.5]GHz, 
b5 = [728.5, 985.6]GHz. Here the background is limited by the 
temperature and emissivity of the cold telescope (total emissivity 
e~ 0.005). 

Experimental Configuration 5 (EC5) is a 4-range differential 
interferometer that operates in deep space (L2), with an actively 
cooled telescope (T < 6K), in the bands bO = [100, 200]GHz, 
bl = [201, 400]GHz, b2 = [401, 600]GHz, b3 = [601, 800]GHz 
with 6 GHz resolution. Here the background is limited by natural 
radiation sources (Galactic, extragalactic, CMB): this allows a 
wider coverage (including high frequencies) without the risk of 
a high background on the detectors. 

For a summary see figH] where the coverage of the different 
configurations is compared to the SZ and foreground spectra. 

3. Simulation of line-of-sight observations 

In order to evaluate the performance of different experimental 
configurations, we have considered the following benchmark sit- 
uation: the observation of a line of sight that crosses a rich cluster 
of galaxies, with t, = 0.005, T e = 5keV, ATcmb = 22fiK (corre- 
sponding to v = 480km/s, if the intrinsic anisotropy of the CMB 
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along the line of sight is negligible), t„, = 0.0001, a = -2.7, 
pi = 2.75x511teV/c). 

We assumed that the angular resolution of all channels of all 
experimental configurations is sufficient to resolve the source. 
Otherwise we would have needed to take into account dilution 
and shape factors: in this LOS approach we avoided these com- 
plications. We believe that this approach is adequate for the pur- 
pose of this paper, which is to compare the performance of dif- 
ferent experimental configurations. 

The integration time on the same line of sight is assumed to 
be 3 hours for all cases but EC4. For EC4 (a space-borne whole- 
sky survey d-la-Planck) the integration time on the considered 
sky pixel is 30 s. 

The signal power on the detectors was computed as 

S(v) = AQ£(v)(l - e,„(v)) [AI t + A/ v + AI CMB + AI nt + AI d ] , (6) 

where A is the collecting area, Q. is the solid angle sampled by 
each detector, E{v) is the efficiency of the detection system, e m (v) 
is the total emissivity of the optical system at room temperature 
and of the atmosphere in the measurement band. 

We have neglected Al/f and AI sy because they are negligible 
with respect to Aid in the frequency range and for the obser- 
vations at high galactic latitudes of interest here. We modelled 
typical Galactic cirrus anisotropy at the angular scale of a clus- 
ter (a few arcmin) as AId{v) = A(v/v ) 4 with A = 600 Jy/sr 
and v = 150 GHz: a value typical of very clean high-latitude 
regions (see e.g. Masi et al. 120061 1. 

We assumed that the detector array is optimized to be limited 
by the photon noise of the radiative background that is produced 
by the instrument and the atmosphere (if present). Cryogenic 
bolometers reach this performance level even in the extremely 
low background achievable in space, if properly designed (see 
e.g. Holmes et al. 120081 . 

The background power on the detectors was computed as 



the FTS the input power is splitted in the two arms of the inter- 
ferometer, and a variable delay is introduced in one of the two 
beams, before recombining them on the detector. The variable 
delay is introduced by a moving mirror, which can be offset by 
x with respect to the corresponding steady mirror in the other 
arm of the interferometer. In this way a 2cx delay is introduced. 
The power measured in position x of the moving mirror is P(x) 
(the interferogram), and the spectrum is estimated as the Fourier 
transform of the interferogram : 



-II 



(P(x) - <P» cos(4jro-x)dx = 



(P(x) - (P)) cos(4no-x)dx, 



(12) 



where <x is the wavenumber (in cm 1 ). The spectral resolution of 
the measurement is 



Act = 1.22/(2x max ) 



(13) 



(see e.g. Chantry [l971l) . In a real instrument the interferogram is 
sampled at the positions xi - iAx, (i - 1..N), with N = x max /Ax. 
The integral is then estimated as a discrete sum 

S (o-j) = 2 Xl^Pixi) - <P» cos(4*o7x,0Ajc = 

= 1.22 jg^PixO - (P))cos(47to-j Xi )/(NAcr). (14) 

Since each sample of the interferogram has an integration time 
T IN, and photons from the whole observed frequency bandwidth 
contribute to the noise, the error in the measurement is therefore 



CTp = 



[!w NEP2 pk dv 



B(v) = AQ£(v) 



e m (v)- 



2hv 3 1 



- 1 



+ [1 - e m (v)] 



2hv i 



1 



- 1 



-(7) 



where x m = hv/kT ln , T m is the temperature of the optical system 
and of the atmosphere, x = hv/kT, T is the temperature of the 
cosmic microwave background. 

The fluctuations of the background were computed as 



y/TTjN 

The error in the measurement of each spectral bin is thus 
a- 2 s = 1 .22 2 Zti o% cos 2 (47ro-x,)/(NAcr) 2 = 

— 1 oo2 °% cos 2 (4;ro-.t,) rv n/\ °% 

~ 1,ZZ NA<t 2 A=l N U -' H /VAo- 2 ' 



or 



(15) 



(16) 



NEP 2 ph =NEPi+NEP 2 CMB 



4k 

NEP 2 m = An—T 5 m E(v)e m (v) 



1 + £(v)e m (v)] 



1) 



(8) 



(9) 



tr s = 0.61— — : = 0.61— -= . (17) 



AcrVr 



cAv Vf 



NEP Z CMB = 

4k 5 , x 4 [e x - 1 + EM(l - e m (v))] 
= Aa^T 5 E(v)(l-e m (y))^ (10) 

For the photometric measurements we integrated equations 
|6]and[8]over the detection bandwidth BW to obtain the error on 
the signal for each band: 



L NEp2 Ph dv 



phot 



2T 



(ID 



where T is the total integration time. 

In the case of spectroscopic measurements, the error on the 
measurement of each spectral bin was computed as follows. In 



For each experimental configuration (EC*) we simulated 
3000 measurements, using equation|6l and adding an error term 
extracted from a Gaussian distribution with zero average and 
standard deviation derived from equations fTTI and [T71 

Typical simulated spectral measurements are reported in fig- 
ure [2] 

4. Results and discussion 

We fitted each simulated measurement using equation [6] In ta- 
ble[T]we report the averages of the best-fit parameters with their 
standard deviation. 

While giving a general idea of the relative efficiency of the 
different configurations, the results reported in table Q] can be 
misleading in the details, since the distributions of the best-fit 
parameters are not Gaussian (nor symmetrical), and there are 
significant correlations between the parameters. This is evident 
from the joint likelihood contours plotted in figures^ 

Earn 
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Fig. 2. From top to bottom: simulated data sets for the spectro- 
scopic configurations EC2 (differential FTS on a stratospheric 
balloon, room temperature), EC3 (differential cold FTS on an 
earth-orbit satellite, with room temperature telescope), EC5 (dif- 
ferential cold FTS on a L2-orbit satellite, with cold telescope). 
The best-fit line through the data points is from eq. [6] The other 
two lines are the thermal plus non-thermal SZ, and the dust 
anisotropy. 



For ECO (4-band ground-based photometer), where only four 
independent data-sets are available for each measurement, we 
tried to fit either three parameters (t,, T, Aid) or four parameters 
(t,, T e , AIj and ATcmb ), adding a fictitious data point with zero 
brightness at zero frequency. The results are dominated by the 
degeneracy between T e and t>, evident from equations Q] and [2] 
without relativistic corrections, the thermal SZ depends on the 



product of electron density and electron temperature. For this 
reason a decrease of, say, a factor 2 of r, is almost perfectly 
compensated for by an increase of a factor 2 of T e . The only 
way to break the degeneracy is through the relativistic correc- 
tions, which, however, are very small: their effect is negligible 
with respect to the typical uncertainties of ground-based mea- 
surements. For this reason we had to use a prior on T e , assuming 
that the information is obtained through independent X-ray mea- 
surements of the specific brightness of the cluster along the same 
line of sight. We tried a very weak prior, with a Gaussian distri- 
bution centred on the true value and a standard deviation of 8 
keV: otherwise the best fit would converge to non-physical val- 
ues for T e . The bias is only mitigated by the introduction of this 
prior. In the 3-parameter case the best fits converge on values 
of the parameters that are not very close (in units of their stan- 
dard deviation) to the input values, and the typical x 1 is high, 
confirming that the effect of the measurement error is negligible 
with respect to the effect of the parameter degeneracies and to 
the necessity of neglecting the non-thermal component in the fit. 

In figure[3]we plot the joint likelihood contours for couples 
of parameters in the 3-parameter fits. The comparison with fig- 
ure [5] shows that part but not all of the bias depends on the 
presence of the non-thermal component: indeed the bias changes 
with Tnt = 0, but is still present. 

If a 4-parameter fit is used, the bias is significantly reduced, 
and does not depend anymore on the presence of the non-thermal 
component. However, this comes at the cost of wider confidence 
areas, as evident from figure [6] 

The contamination of dust is significant and degenerates to 
some extent also with ATcmb- in the 3-parameter fit, where 
ATcmb is not fitted, the best-fit for Aid is far from the input value, 
and this probably contributes to the tension between the best-fit 
value of T e and its input value. We believe that this is the result 
of a degeneracy of parameters combined with the poor cover- 
age of high frequencies in ground-based experiments. In the 4- 
parameter fits the x 2 improves, but the best fits are still biased, 
and only with the 3 keV prior the estimates of all parameters are 
consistent with the input values (see fig.|6j- 

For EC1 (4-bands balloon-borne photometer) the coverage 
of frequencies higher than 300 GHz helps in removing the de- 
generacy (see fig. |7). If we use the very weak prior on T e 
(Gaussian with standard deviation 8 ke V), the best-fit cluster pa- 
rameters are already close to the input values, even if the error 
(completely dominated by the degeneracy) is relatively large. 
The situation improves, of course, if the standard deviation on 
the prior is reduced to 3 keV. ATcmb is better constrained, but is 
biased low. While performing significantly better than ECO, this 
configuration is limited by radiation noise in the high-frequency 
bands, which can be removed only by cooling the telescope 
or/and adding spectroscopic capabilities (see below). 

In fig. [8] we plot the histograms of the best-fit T e x t> nor- 
malized to the input of the simulation. It is evident how the bias 
is reduced passing from ECO with three parameters to ECO with 
four parameters to EC1, which shows how important the cover- 
age of high frequencies is, which are difficult to observe from 
the ground. 

For EC4, where six independent data-sets are available for 
each measurement, we tried to fit four parameters (t>, T e , AI^ 
and ATcmb, see fig. [9]) and also six parameters (including t„, 
and p\, see fig. [Toll, always adding the fictitious zero frequency 
zero brightness data point. The good coverage of high frequen- 
cies and low photon noise owing to the low radiative background 
results in very good performance in terms of statistical errors. 
Again, the presence of the non-thermal component produces 
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0.000 I , 

-0.003 



A I. (MJy/sr) T.(keV) 

Fig. 3. Joint likelihood contours (95.6% and 68.7%) for couples 
of best-fit parameters for the experimental configuration ECO 
(ground based observations with a 4-band photometer) where the 
data are analysed by fitting the three parameters r t , T e , Aid- The 
filled contours are for a Gaussian prior on T e with cr = 3 keV; the 
dashed contours are for a Gaussian prior with cr — 8 keV. The * 
symbols mark the input values of parameters. In this case a sig- 
nificant bias for all parameters is evident. The anti-correlation 
between r, and T e is also evident. 



Fig. 6. Same as fig. [3] when the data are analysed by fitting four 
parameters (r,, T e , Aid and ATcmb)- The bias is reduced, but 
the anti-correlation between t, and T e is still evident and the 
confidence contours are wider. 






-0.003 0.000 0.003 0.1 

A I, (MJy/sr) 

Fig. 4. Same as fig. [3] but with t„ 
does not vanish. 



12.5 25.0 
T.(keV) 

= 0. The bias changes, but 



Fig. 7. Joint likelihood contours (95.6% and 68.7%) for couples 
of best-fit parameters for the experimental configuration EC1 
(balloon-borne observations with a 4-band photometer, extend- 
ing the coverage to high frequencies not observable from the 
ground) where the data are analysed by fitting four parameters 
(t,, T e , Aid, ATcmb)- The filled contours are for a Gaussian prior 
on T e with cr — 3 keV; the dashed contours are for a Gaussian 
prior with cr = 8 keV. The * symbols mark the input values of 
parameters. 



6 



de Bernardis et al. : Cluster parameters with SZ 




1.0 1.1 1.2 1.3 

(T„ x t, ) m / (T e x r, ) lnput 

Fig. 8. Histograms of the best-fit T e x t, normalized to the in- 
put of the simulation, for ECO with three parameter fits (dashed 
lines), ECO with four parameter fits (continuous lines) and EC1 
with four parameter fit (dot-dashed lines). Thin lines are for a 
Gaussian prior on T e centred on the input value of T e and with a 
standard deviation of 8 keV; thick lines are for a standard devia- 
tion of 3 keV. 



a bias in the determination of the other parameters for the 4- 
parameter fit. In the 6-parameter fits, the non thermal parameters 
are basically not constrained (with a bimodal distribution of the 
best fit p\), but the other parameters are well constrained and un- 
biased. The effects of parameters degeneracies are still evident, 
however, and it is difficult to estimate ATcmb because it tends to 
be biased low in the 4-parameter fit, while is not well constrained 
in the 6-parameter fits (see also tableQ]i. 

The space spectrometers EC2, EC3, and EC5, featuring wide 
frequency coverage and low radiative background noise, perform 
better, and allow an unbiased recovery of six parameters with in- 
creasing accuracy. If we compare the parameter space volume 
constrained by the different configurations, we find that with re- 
spect to EC2, EC3 reduces the volume by a factor 4.6, and EC5 
by a factor 27 (with the 8 keV prior). 

The balloon-borne warm spectrometer EC2 is still limited 
at high frequencies by radiative background fluctuations. For 
this reason ATcmb is basically not constrained (see fig. ITTb . 
However, the other parameters are unbiased and well con- 
strained even with the very weak prior on T e . An experiment 
like OLIMPO, which combines photometric measurements (as 
in EC1) and spectroscopic measurements (as EC2), performing 
both during the same flight, can use the first measurement to 
optimally constrain Aid and ATcmb, and the second to better 
constrain r, and T e . Note that even r nt is close to be detected 
(and can be detected with an integration time longer than the 3 
hours considered here) and p\ is also constrained (still with a 
bimodal distribution). Finally note that the spectroscopic capa- 
bilities allow the user to remove the contamination of Galactic 
CO and other cooling lines within the photometric bands unam- 
biguously. 

The low-Earth-orbit case EC3 improves over the balloon 
spectrometer EC2. Cooling the spectrometer results in a signif- 
icant unbiased measurement of all six parameters, including the 
elusive ATcmb (see fig.fT2l). Because the noise is reduced, the ef- 
fect of degeneracies is more evident. The 3 keV prior is needed 
to fully exploit the potential of this configuration. 




Fig. 9. Joint likelihood contours (95.6% and 68.7%) for couples 
of best-fit parameters for the experimental configuration EC4 (6- 
band cryogenic photometer operating in L2) where the data are 
analysed by fitting four parameters (t>, T e , Aid, ATcmb)- The 
filled contours are for a Gaussian prior on T e with cr — 3 keV; 
the dashed contours are for a Gaussian prior with cr = 8 keV. 
The * symbols mark the input values of parameters. 



Experimental configuration EC5, where both the spectrome- 
ter and the telescope are cold, results in an improved precision 
of the determination of all parameters, and little sensitivity to the 
prior on T e (see fig.[T3land table[TJ. 

We have investigated the possibility of extending the fre- 
quency coverage of EC5 up to 1 THz. Keeping the same simpli- 
fied model for dust, we found an improvement of about a factor 2 
in the uncertainties of all parameters but ATcmb (which does not 
improve) and Aid (which improves by a factor 5). However, our 
simplistic model with a single dust population probably starts 
to be insufficient at these high frequencies. There are, however, 
other drivers to extend the frequency coverage of such an am- 
bitious experiment as high as possible: the study of interstellar 
atoms and molecules and the study of galaxies in the redshift 
desert would benefit significantly from such an extension. 



5. Conclusions 

With this simulation-based study we confirmed that few-band 
photometric measurements of the SZ effect can be significantly 
biased, and that most of the bias is removed by adding bands and 
covering high frequencies (v £ 240GHz). This requires the use 
of space-borne experiments. Low-resolution spectroscopic mea- 
surements of the SZ are extremely promising. An exploratory 
experiment, featuring a four-band photometer on a balloon and 
a plug-in room-temperature FTS [configurations EC1 and EC2, 
modelled on the OLIMPO experiment (Masi et al. I2008H is re- 
quired to confirm that the expectations estimated here can really 
be achieved. In other words, we relied here on the differential 
properties of FTS instruments: we need to demonstrate experi- 
mentally that systematic effects are negligible even with the ex- 
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f 10 
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Fig. 10. Same as fig. [9] (configuration EC4) when the data are 
analysed by fitting six parameters (t,, T e , Aid, &T C mb, T nt , pi). 



Fig. 12. Same as fig.[TT|for the experimental configuration EC3 
(cold spectrometer on an Earth-orbit satellite with radiatively 
cooled (80K) telescope) 



Fig. 11. Joint likelihood contours (95.6% and 68.7%) for cou- 
ples of best-fit parameters for the experimental configuration 
EC2 (4-range warm spectrometer on a balloon platform). The 
filled contours are for a Gaussian prior on T e with cr — 3 keV; 
the dashed contours are for a Gaussian prior with cr = 8 keV. 
The * symbols mark the input values of parameters. 



treme requirements of SZ measurements. OLIMPO, in the dual 
photometric-spectroscopic configuration, is a perfect demonstra- 
tor, and is expected to produce exciting improvements on the 
physics of cosmic structures derived from the SZ effect. Satellite 
missions can accommodate a cryogenic FTS (a large cryostat 



Fig. 13. Same as fig.[TT|for the experimental configuration EC5 
(cold spectrometer on a satellite in L2, with actively cooled (4K) 
telescope). 



is required anyway to provide the required hold time), reducing 
the instrumental background at high frequency. A Molniya-orbit 
mission [like the proposed SAGACE (de Bernardis et al. 2010)] 
represents a cost-effective solution, able to provide an exten- 
sive catalogue of unbiased SZ measurements in various cosmic 
structures (e.g. galaxy clusters and groups, radiogalaxy lobes, 
gaseous halos of galaxies). A cold telescope in L2 (like the pro- 
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posed Spectrum-M - millimetron mission) would be limited only 
by the intrinsic degeneracy of the parameters, and would open 
new horizons in SZ science, allowing a precise determination of 
both the thermal and non-thermal components of the plasma in 
the cluster. 
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Parameter 


input 


best fit ECO 


best fit ECO 


best fit ECO 


best fit ECO 


best fit EC1 


best fit EC1 






ground - phot. 


ground - phot. 


ground - phot. 


ground - phot. 


balloon - phot. 


balloon - phot. 






prior cr= 3 keV 


prior cr= 8 keV 


prior cr= 3 keV 


prior cr= 8 keV 


prior <r= 3 keV 


prior <x= 8 keV 


r, x 10 j 


5 


5.2 ± 0.6 


3.9 ± 1.3 


5.1 ±0.7 


4.7 ±3.1 


5.0 ± 1.3 


5.3 ± 1.7 


T(keV) 


8.5 


10.8 ± 1.2 


15.8 + 3.5 


8.7 ± 1.1 


12±4 


8.5 ±0.3 


8.6 ± 1.7 


AT cmb (mK) 


22 






19 ± 6 


11 ±6 


16 ±2 


16 ±4 


hWylsr) 


600 


-1800 ±200 


-1100 ±500 


390 ± 480 


50 ± 660 


580 ± 70 


590 ± 80 


T„, X 10 J 


0.1 














pi(511 keV/c) 


2.75 














(X 2 )/D0F 




6.3/2 


4.3/2 


1.0/1 


1.0/1 


1.1/1 


1.1/1 



Parameter 


input 


best fit EC4 


best fit EC4 


best fit EC4 


best fit EC4 






L2 - cold phot. 


L2 - cold phot. 


L2 - cold phot. 


L2 - cold phot. 






prior cr= 3 keV 


prior cr= 8 keV 


prior cr= 3 keV 


prior <x= 8 keV 


T, x 10 J 


5 


5.2 ± 0.9 


5.4 ± 1.9 


5.0 ±0.8 


5.3 ±2.1 


T(keV) 


8.5 


8.5 ± 1.5 


8.7 ±2.5 


8.6 ± 1.3 


9.0 ± 2.7 




22 


16±4 


16±4 


21 ±22 


22 ±23 


AWylsr) 


600 


601 ±6 


600 ±6 


601 ±5 


601 ±6 


T,„ x 10 J 


0.1 






0.1 ±0.4 


0.1 ±0.4 


pi(511 keV/c) 


2.75 






6±4 


6±4 


( X z )/DOF 




2.7/3 


2.4/3 


1.3/1 


1.0/1 



Parameter 


input 


best fit EC2 


best fit EC2 


best fit EC3 


best fit EC3 


best fit EC 5 


best fit EC 5 






balloon - warm spec. 


balloon - warm spec. 


EO - cold spec. 


EO - cold spec. 


L2 - cold spec. 


L2 - cold spec. 






prior cr= 8 keV 


prior cr= 3 keV 


prior <x= 8 keV 


prior <r= 3 keV 


prior <r= 8 keV 


prior <r= 3 keV 


T, x 10 j 


5 


5.0 ±0.9 


4.9 ±0.8 


5.8 ±2.6 


5.2 ± 0.6 


5.1 ±0.6 


5.1 ±0.5 


T(keV) 


8.5 


8.4 ±0.8 


8.5 ±0.1 


7.7 ± 2.0 


8.1 ±0.8 


8.5 ± 1.2 


8.5 ± 1.0 


ATcmb(mK) 


22 


20 ± 50 


20 ± 50 


23 ±8 


22 ± 8 


22 ±4 


22 ±4 


AWylsr) 


600 


570 ± 270 


560 ± 270 


590 ± 40 


590 ± 40 


600 ±4 


600 ±4 


T„, x W 


0.1 


0.1 ±0.1 


0.1 ±0.1 


0.12 ±0.03 


0.11 ± 0.02 


0.10 ± 0.01 


0.10± 0.01 


pi(511 keV/c) 


2.75 


2.6 ±0.7 


2.5 ±0.7 


2.5 ±0.9 


2.7 ± 1.1 


3.0 ± 1.0 


2.9 ±0.9 


(X^IDOF 




34.9/34 


34.9/34 


77.8/78 


78.0/78 


110.0/110 


110.1/110 



Table 1. Results of simulations. For each experimental configuration (EC*) we simulated 3000 measurements and fitted each 
simulated measurement using equation|6] We report the averages of the best-fit parameters with their standard deviation. In the top 
table we summarize the results for 4-band photometers for both weak and medium prior for T e ; in the middle table we report the 
results for a 6-band photometer (Planck-HFI bands). In the bottom table we report the results for low-resolution (Av = 6 GHz) 
spectrometers with increasing mission complexity (balloon, Earth orbit, orbit around the Sun-Earth Lagrangian point L2). 
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